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2000-2001  Annual  Progress  Report  (Year  1) 


This  report  presents  the  specific  aims  and  accomplishments  of  our  breast  cancer  research 
project  during  the  first  year  of  funding  sponsored  by  the  U.S.  Department  of  the  Army.  It  covers 
our  activities  from  September  1, 2000  to  October  31, 2001. 

Introduction 

The  goal  of  this  research  project  is  to  develop  and  evaluate  a  new  approach  to  monitoring  of 
vascular  blood  oxygen  saturation  (SO2)  and  blood  concentration  (Hb)  of  breast  tumors  using  near 
infrared  (NIR)  spectroscopy  and  imaging  techniques.  Further,  we  wish  to  compare  and  validate  the 
optical  method  with  concurrent  measurements  of  tumor  oxygen  tension  on  the  basis  of  19F  EPI 
relaxation  of  hexafluorobenzene  (HFB).  Once  we  have  rigorously  established  the  integrated 
approach  to  tumor  physiology,  we  will  explore  the  influence  of  three  interventions  expected  to 
modify  tumor  physiology.  A  better  understanding  of  the  interplay  of  these  parameters  in  the  natural 
history  of  a  tumor  could  lead  to  enhanced  therapeutic  approaches  and  provide  a  novel 
diagnostic/prognostic  tool  for  breast  cancer  research  and  clinical  practice. 

The  project  has  four  specific  aims: 

Aim  1:  To  evaluate  a  single-channel,  dual  wavelength,  NIR,  frequency-domain  oximeter  and  the 
algorithms  for  obtaining  tumor  SO2  against  tumor  p02  measured  by  l9F  EPI  relaxation  of  HFB. 

Aim  2:  To  modify  the  single-channel  system  into  a  3-channel  NIR  system. 

Aim  3:  To  investigate  heterogeneity  of  SO2  in  breast  tumors  using  the  3-channel  NIR  system. 

Aim  4:  To  study  the  influence  of  three  interventions  on  SO2  and  pC>2  of  the  tumor  simultaneously. 

Specifically,  Task  1  was  planned  for  months  1-15  to  accomplish  Aim  1 : 

Task  1:  To  evaluate  a  single-channel,  dual  wavelength,  NIR,  frequency-domain  oximeter  and  the 
corresponding  algorithm  for  obtaining  tumor  SO2  against  tumor  p(>2  measurement  using  l9F  MR 
EPI  of  hexaflourobenzene,  (months  1-15): 

a.  test  and  calibrate  the  existing  single-channel,  dual  wavelength,  NIR  oximeter  using 
blood-yeast  models  to  achieve  a  stable  and  low-noise  system; 

b.  perform  laboratory  NIR  and  MR  EPI  measurements  on  rat  breast  tumors  under  two 
respiratory  conditions  (33%  oxygen  as  the  base  line  and  100%  as  an  intervention). 

c.  employ  proper  algorithms  to  calculate  absolute  values  and  changes  of  SO2  and  Hb  in 
the  tumor  vascular  bed  before,  during,  and  after  the  respiratory  challenge. 

d.  Compare  the  SO2  values  with  pC>2  and  study  vascular  and  tissue  dynamic  response  to 
the  respiratory  change. 


Body  of  the  Report 

As  mentioned  in  Introduction  section,  our  Year  1  focus  was  to  achieve  Task  1,  and  we  have 
mainly  accomplished  this  task,  as  reported  below: 

Task  la:  to  test  and  calibrate  the  existing  single-channel,  dual  wavelength.  NIR  oximeter 
using  blood-yeast  models  to  achieve  a  stable  and  low-noise  system: 

The  PI  and  her  graduate  student,  Jae  Kim,  have  made  significant  efforts  to  accomplish  this 
task.  Two  kinds  of  experiments  were  performed  to  calibrate  our  NIR  system:  1)  cross-talk 


measurement  between  amplitude  and  phase  of  the  instrument  and  2)  blood  phantom 
measurement  for  error  analysis. 


Crosstalk  Study 

In  principle,  changes  in  amplitude  of  a  frequency-domain  photon  migration  system  are 
independent  of  its  changes  in  phase,  i.e.,  the  crosstalk  between  amplitude  and  phase  should  be 
minimal  for  an  ideal  NIR  system.  However,  in  reality,  it  is  difficult  to  have  such  a  perfect 
condition  for  an  actual  instrument  and  necessary  to  calibrate  the  system.  To  conduct  the 
crosstalk  test,  we  altered  optical  filters  in  front  of  the  two  light  sources  (758  nm  and  785  nm)  to 
attenuate  the  optical  amplitude,  so  we  could  see  whether  phase  data  were  influenced  by 
amplitude  changes.  The  test  showed  that  the  phase-amplitude  crosstalk  exists  if  the  measured 
output  signal  is  either  too  large  or  too  small.  We  concluded  that  the  output  range  of  our 
instruments  should  be  between  350  mV  and  100  mV  during  measurement  (Figs.  1  and  2)  to 

avoid  the  unwanted  phase-amplitude  crosstalk. 
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Figure.  1  Crosstalk  test  for  75 8nm  wavelength 


Figure.  2  Crosstalk  test  for  785nm  wavelength 


We  also  plotted  O.D  (optical  density)  of  the  filters  versus  log(amplitude  of  the  measured 
output  signal)  to  see  the  electrical  response  of  the  instrument  when  the  input  light  intensity 
changes.  It  was  found  that  as  long  as  we  are  in  the  output  range  of  100  mV  to  350  mV,  we  can 
achieve  good  linearity  for  the  system  response,  as  shown  in  Figures  3  and  4. 


O.D.  O.D. 


Figure  3.  Linear  system  response  between  O.D.  and  Figure.  4  Linear  system  response  between  O.D.  and 

Log(amplitude)  (starts  from  350mV  amplitude).  Log(amplitude)  (starts  from  193mV  amplitude). 


Blood  Phantom  Study 


To  calibrate  our  NIR  system  for  quantifying  hemoglobin  concentrations,  we  used  human 
blood  for  tissue  phantom  experiments.  The  measurements  were  conducted  at  the  University  of 
Pennsylvania  since  the  graduate  student,  Jae  Kim,  needed  some  training  in  conducting  this  kind 
of  experiments.  One  packet  of  SIGMA  P-3813  phosphate  buffered  saline  (pH  7.4)  powder  was 
used  to  make  a  2-litter  buffer,  and  100  ml  of  20%  intralipid  were  added  into  the  buffer  for  a  1% 
intralipid  solution.  Then  14g  of  baking  yeast  were  mixed  with  tap  water  before  being  added  into 
the  2-litter  buffer  solution. 

Amount  of  blood  for  each  addition  to  the  solution  was  3  ml  for  the  first  two  additions  and  6 
ml  for  later  additions.  For  the  first  two  additions,  we  added  3  ml  of  blood  each  time,  which 
contained  6.6  pmoles  of  hemoglobin  and  gave  a  concentration  of  3.3  pM  in  the  solution  (=  6.6 
pmoles  /  2  L).  The  later  6  ml  of  blood  addition  resulted  in  an  addition  of  6.6  pM  (=  13.2 
pmoles/2L)  in  hemoglobin  concentration. 

In  the  NIR  range,  major  light  chromophores  in  tissues  are  Hb  and  Hb02.  Therefore,  the 
absorption  coefficients  at  two  wavelengths  can  be  associated  with  the  concentrations  of  Hb  and 
Hb02,  expressed  as  [Hb]  and  [Hb02],  respectively,  by 

pa758  =  £Hb758[Hb]  +  £Hb02758[Hb02],  (1) 

Pa785  =  SHb785[Hb]  +  SHb02785[Hb02]  .  (2) 

Our  algorithm  is  based  on  modified  Beer-Lambert’s  Law,  which  has  the  following  expression. 

ln(I0/I)  =  =  paL,  (3) 

where  I0  and  I  are  the  incident  and  the  detected  optical  intensities,  respectively,  and  L  is  the 
optical  path  length  traveled  by  light  inside  the  medium. 

Therefore,  changes  of  absorption  coefficients  can  be  expressed  as 


Apa  paT  ■  paB  ln(  I0/I|)L  -  ln(I0/IB)L  In  (Ig/Ij)  /  L,  (4) 

with  the  assumption  that  L  is  constant.  Then  the  changes  of  Hb  and  Hb02  concentrations  can  be 
calculated  by 

7C0  705 

A[Hb02]  =  [-5.73*ln(IB/IT)  +  9.67*ln(IB/IT)  ]/L,  (5) 

758  785 

A[Hb]total  =  A[Hb]  +  A[Hb02]=[-l  ,32*ln(IB/IT)  +  6.23*ln(IB/IT)  ]  /  L,  (6) 

where  Ib  =  baseline  amplitude;  IT  =  transition  amplitude;  L  =  optical  pathlength  between 
source/detector;  the  constants  were  computed  with  the  extinction  coefficients  for  human  oxy  and 
deoxyhemoglobin  at  the  two  wavelengths  used.  Some  details  for  the  algorithms  can  be  found  in 
the  manuscript  attached  in  Appendix  a. 
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By  using  eqs.  (5)  and  (6),  we  can  obtain  [Hb02]  and  [Hb]totai,  as  plotted  in  Figure  5.  This 
figure  basically  shows  that  33  pM  and  30  pM  of  [Hb]totai  increased  for  the  first  two  additions  of  3 
ml  blood,  and  52  pM  and  59  pM  of  [Hb]totai  increased  for  the  third  and  fourth  6  ml  blood 
additions.  These  are  about  10  and  8.5  times  larger  than  3.3  pM  and  6.6  pM,  respectively. 
Further  studies  are  being  conducted  to  understand  the  origin  of  this  large  factor. 


Yeast  02  Yeast  02  Yeast  02  Yeast 


Figure  5.  Phantom  experiment  using  1%  intralipid  in  a  saline  solution.  Changes  of  [Hb],  [Hb02],  [Hb]totai  during 
oxygenation  and  deoxygenation  process  were  measured  by  using  N1RS. 


During  the  studies,  we  also  learned  that  a  small  variation  in  extinction  coefficients  could 
cause  large  variations  in  determinations  of  [Hb02]  and  [Hb]totai-  Extinction  coefficients  of  human 
oxy-hemoglobin  and  deoxy-hemoglobin  can  be  found  in  two  sources  given  by  Zijlstra  and 
Cope1,2.  Until  now  we  had  used  Zijlstra’s  extinction  coefficients,  but  we  have  found  that  using 
Mark  Cope’s  extinction  coefficients  give  us  better  results  of  changes  in  hemoglobin  oxygenation. 
We  applied  Cope’s  extinction  coefficients  in  this  report.  In  Mark  Cope’s  dissertation,  extinction 
coefficients  of  hemoglobin  were  quoted  for  a  molecular  weight  of  64450  and  hence  four  times 
larger  than  the  values  that  are  quoted  as  “per  equivalent”.  Each  hemoglobin  molecule  consists  of 
four  protein  chains  called  globins,  each  of  which  is  bound  to  one  heme,  a  red-pigmented 
molecule  that  contains  iron.  Therefore,  to  use  his  extinction  coefficients  as  “per  equivalent”  we 
need  to  divide  his  extinction  coefficient  by  four. 

The  changes  in  total  hemoglobin  concentration  in  the  blood  phantom  study  showed  that  Mark 
Cope’s  extinction  coefficients  divided  by  4  resulted  in  the  theoretical/expected  values  when  we 
assumed  that  the  DPF  of  the  phantom  solution  is  about  9  (please  refer  to  the  manuscript  in 
Appendix  a  for  DPF).  We  are  continuing  to  explore  the  crucial  relationship  between  the 
extinction  coefficients  and  accurate  determinations  for  [Hb02]  and  [Hb]totai* 
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Task  lb:  to  perform  laboratory  NIR  and  MR  EPI  measurements  on  rat  breast  tumors  under 
two  respiratory  conditions  (33%  or  21%  oxygen  as  the  base  line  and  100%  or  Carbogen  as  an 
intervention); 

Basically,  the  manuscript  attached  in  Appendix  a)  describes  both  the  methodology  and 
results  for  this  sub-task.  Please  refer  to  it  for  detailed  information.  Note  that  although  the  study 
reported  in  the  manuscript  used  primarily  prostate  tumors,  the  methodology  and  the  algorithm 
developed  in  the  study  are  applicable  to  breast  tumors.  The  actual  ivestigations  for  breast  tumors 
are  underway. 


Task  lc:  to  employ  proper  algorithms  to  calculate  absolute  values  and  changes  of  SO?  and 
Hb  in  the  tumor  vascular  bed  before,  during,  and  after  the  respiratory  challenge; 

The  PI  and  her  graduate  student,  Mengna  Xia,  have  been  working  on  this  sub-task,  and  the 
status  report  for  this  is  as  follows: 

In  order  to  quantify  Hb  and  SO2  of  a  measured  sample,  we  have  attempted  two  algorithms. 


Algorithm  1  (variation  method):  According  to  Matthias  Kohl  et  al.3’4,  small  changes  in 
absorption,  pa,  induce  changes  in  both  measured  reflectance,  R,  and  phase,  <1>,  and  they  can  be 
expressed  as 


-Pap 


R(P)  =  za& /P  +  Mcff)- — T 
l7tp 


<t>(p)  = 


2  c(D  +  pJpj5) 


(7) 


<f>  =  -2nv  <  t  > 


where  p  is  the  source-detector  separation,  zo  =l/(pa  +  Bs’X  M-eff=[3pa(Pa+Ps’)]°  5>  D=T/(3ps’),  pa 
and  [is  are  absorption  and  reduced  scattering  coefficients,  respectively,  t  is  time,  and  v  is  the 
modulation  frequency  of  the  diffuse  photon  density  waves  traveling  in  tissues.  Based  on  eq.  (7) 
and  defining  A  as  A=log(R),  we  derived  an  equation  to  obtain  pa: 


Pc, 


1  + 


12  p2p:qc, 


a 


6  p2p.: 


(8) 


where  a=c/[7t:vln(10)]  and  Qa  =  AA /AO .  Qa  can  be  obtained  by  measuring  changes  in  amplitude, 

AA,  and  phase,  AO.  The  condition  to  employ  eq.  (8)  to  quantify  pa  is  that  the  measured  sample 
has  to  undergo  a  perturbation  so  that  small  changes  in  both  AA  and  AO  occur  and  lead  to  Qa 
determination.  Further  inspection  on  Qa  shows  that  Qa  is  largely  independent  of  the  scattering 
property,  ps\  of  the  medium,  as  shown  in  Figure  6  below.  This  means  pa  obtained  from  Qa  is 
largely  independent  of  ps’  variations. 
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Qa  vs.  Ua  (Us'=5,10,15(l/cm)) 


Figure  6.  Qa  plot  versus  pa  at  three 
different  ps’  values  of  5,  10,  and  15  cm'1. 
The  differences  of  Qa  values  with 
different  ps’  are  minimal. 


Algorithm  2  (Iterative  fitting):  On  the  other  hand,  based  on  the  diffusion  approximation  with  a 
semi-infinite  boundary  condition5,  we  developed  an  iterative,  curve-fitting  program  for  our 
single  channel  (single  source-detector  separation),  two  wavelength,  IQ  system  to  calculate  pa  and 
ps’  values  of  the  measured  samples.  The  prerequisite  for  this  algorithm  to  be  applied  is  that  we 
have  to  know  one  pair  of  pa  and  ps’  for  calibrating  the  instrument.  In  the  following  is  the  protocol 
for  the  curve-fitting  algorithm  in  getting  the  pa  and  ps’  values  at  different  times  (or  under 
different  conditions): 

a)  to  have  a  set  of  known  pao  and  ps0’  for  calibration; 

b)  to  have  an  initial  guess  for  a  set  of  pa  and  ps’  to  calculate  the  reflectance,  R(p,o),pa,ps’); 

c)  to  compare  the  calculated  R(p,to,  pa,  and  ps’)  with  the  real  measured  R,  and  try  to 
decrease  the  difference  between  two  R’s  by  iteratively  changing  the  values  of  pa  and  fV- 
The  set  of  pa  and  ps’  that  give  the  minimal  difference  between  the  two  R’s  is  the  best  fit 
to  the  diffusion  equation,  and  thus  it  is  considered  to  be  the  pa  and  ps’  values  of  the 
measured  sample. 


Steps  a)  to  c)  can  be  repeated  for  different  wavelengths,  and  the  determinations  of  pa  at 
two  wavelengths  result  in  the  calculations  of  absolute  concentrations  of  [Hb],  [Hb02],  and 


[Hbtotai],  which  are  derived  and  expressed  as  follows: 

[Hb]=1.1035  pa758  -  0.8  5  83pa785  (9) 

[Hb02]=  -1.4331  pa758  +  2.4165  pa785  (10) 

S02=[Hb02]/(  [Hb02]+[Hb])  (11) 

[HbJtota.  =  [Hb02]+[Hb]  (12) 


The  following  figures,  Figures  7-9,  show  an  example  of  quantifying  pa,  ps’,  and  S02  from  the 
same  blood  phantom  experiment  whose  result  was  shown  above  in  Figure  5: 


Figure  7.  Absorption  coefficient  plot 
versus  time  at  two  wavelengths,  758nm 
and  785  nm. 
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Figure  8.  Reduced  scattering  coefficient 
plot  versus  time  at  two  wavelengths, 
758nm  and  785  nm. 
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For  our  breast  tumor  studies,  our  instrument  has  only  one  source-detector  separation  and 
two  wavelengths,  so  we  have  to  initially  assume  a  ps’  value  in  order  to  obtain  our  initial  pa  value, 
basic  on  Algorithm  1 .  Furthermore,  we  know  that  ps’  values  of  biological  tissues  are  in  the 
range  from  5  cm-1  to  15  cm-1.  For  such  a  range,  through  computer  simulations,  we  learned  that 
both  Algorithm  1  and  Algorithm  2  mentioned  above  are  not  too  sensitive  to  variations  of 
scattering  properties  when  calculating  total  hemoglobin  concentrations.  However,  both 
algorithms  have  their  limitations.  For  Algorithm  1,  it  can  give  only  the  pa  value  in  the  period 
with  a  perturbation  in  pa,  which  is  induced  by  changes  in  hemoglobin  concentration.  For 
Algorithm  2,  it  can  give  pa  values  only  after  one  set  of  pa  and  ps’  is  known  since  we  have  to  use 
them  for  calibrating  our  instrument. 

A  Combined  Algorithm: 

Thus,  we  combined  these  two  algorithms  to  obtain  pa  of  the  measured  sample  during  a  period 
of  intervention/perturbation.  In  the  following  are  listed  two  basic  steps: 

a)  to  use  Algorithm  1  (i.e.,  the  variation  method)  to  get  one  pa  value  within  a  perturbation 
period,  and  assume  ps’=T0  cm'1; 

b)  then,  to  use  this  determined  pa,  along  with  ps’=10  cm'1,  as  the  calibration  value,  and  apply 
them  in  Algorithm  2  (i.e,  the  iterative  fitting)  to  get  pa  values  of  the  other  points  in  the 
measured  period. 

Our  studies  through  computer  simulations  have  shown  that  the  errors  caused  by  the 
assumption  of  ps’=10  cm'1  are  within  10-15%  if  the  true  ps’  values  of  the  measured  samples  are 
in  the  range  of  5-15  cm'1  Right  now,  we  are  conducting  laboratory  experiments  to  validate  the 
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combined  algorithm  for  Hb  and  Hbtotai  calculations  and  having  a  manuscript  written  in  1-2  month 
to  describe  the  combined  algorithm. 

Task  Id:  to  compare  the  SO?  values  with  pO?  and  study  vascular  and  tissue  dynamic  response 
to  the  respiratory  change. 

Again,  the  manuscript  attached  in  Appendix  a)  describes  both  the  methodology  and  results 
for  this  sub-task.  Please  refer  to  it  for  detailed  information.  Note  that  although  the  tumors  used 
in  the  study  were  primarily  prostate  tumors,  the  methodology  and  the  algorithm  developed  in  the 
study  are  applicable  to  breast  tumors.  The  actual  studies  for  breast  tumors  are  underway. 


Key  Research  Accomplishments 

1)  We  have  investigated  the  relationships  between  tumor  vascular  oxygenation  and  tumor 
oxygen  tensions  under  carbogen  intervention  using  the  NIR  spectrometer,  p02  needle 
electrodes,  and  19F  MRI  p02  mapping.  This  study  demonstrates  that  the  NIRS 
technology  can  provide  an  efficient,  real-time,  non-invasive  approach  to  monitoring 
tumor  physiology,  and  emphasizes  the  need  to  develop  an  imaging  technique  to  reveal 
spatial  heterogeneity. 

This  part  of  work  has  been  mainly  carried  by  Jae  Kim,  a  Ph.D.  graduate  student  in  our 
BME  program. 

2)  We  have  developed  a  combined  algorithm  that  can  be  used  to  quantify  oxygenated  and 
total  hemoglobin  concentrations  based  on  our  single  source-detector  separation,  two- 
wavelength,  frequency-domain  NIR  system.  This  algorithm  will  be  applied  to  breast 
tumor  measurements  for  absolute  quantifications. 

This  part  of  work  has  been  mainly  carried  by  Mengna  Xia,  a  Ph.D.  graduate  student  in 
our  BME  program. 

3)  We  have  performed  a  few  preliminary  studies  on 

a.  determination  of  breast  tumor  oxygen  consumption  by  NIRS  , 

b.  vascular  oxygen  dynamics  of  breast  tumors  in  response  to  pharmacological 
interventions  by  NIRS, 

c.  blood  volume  quantification  in  breast  tumors  by  19F  MRS  of  PFOB  and  its 
correlation  to  A[Hb]T0tai  by  NIRS. 

This  part  of  work  was  mainly  conducted  by  Yulin  Song ,  who  was  a  Ph.D.  candidate  in 
our  Biomedical  Engineering  Program  from  1995  to  2001.  He  completed  his  Ph.D.  work 
in  August  of  2001  with  a  dissertation  title  of  “Tumor  Oxygenation  and  Physiology  in 
Response  to  Therapeutic  Interventions  Investigated  by  19 F  MRI  an  NIR  Spectroscopy.  ” 
See  appendix  b  for  details. 


Reportable  Outcomes 

Manuscripts  for  Peer-reviewed  Journals: 

1.  Jae  G.  Kim,  Yulin  Song,  Dawen  Zhao,  Anca  Constantinescu,  Ralph  P.  Mason,  and  Hanli 
Liu,  “Interplay  of  Tumor  Vascular  Oxygenation  and  Tumor  p02  Observed  Using  NIRS, 
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Oxygen  Needle  Electrode,  and  19F  MR  p02  Mapping,”  submitted  to  J.  of  Biomedical 
Optics  (2001). 

2.  Jae  Kim  and  Hanli  Liu,  “Error  Analysis  of  Calculated  Elemoglobin  Concentrations 
Caused  by  Variations  of  Extinction  Coefficients,”  in  preparation  and  to  be  submitted  to 
Optics  Express. 

3.  Mengna  Xia  and  Hanli  Liu,  “Determination  of  Total  Hemoglobin  Concentrations  in 
Blood-Perfused  Turbid  Meida  with  Scattering  Variations,”  in  preparation  and  to  be 
submitted  to  J.  of  Biomedical  Qptjcs- 

Presentations  and  Proceeding  Papers: 

1 .  Hanli  Liu,  “Near-Infrared  Spectroscopy  Used  in  vivo  for  Assisting  Neurosurgery  and 
Tumor  Oximetry,”  Workshop  on  Biotechnology  and  Nano-biotechnology — Challenges 
and  Opportunities  of  Drug  Targeting  in  the  Post  Genomic  Era,  organized  by  the 
University  of  Texas  at  Dallas  and  the  University  of  Texas  Southwestern  Medical  Center 
at  Dallas,  Dallas,  October  18,  2001. 

2.  J.  G.  Kim,  Y.  Song,  D.  Zhao,  A.  Constantinescu,  R.  P.  Mason,  and  Hanli  Liu,  “Interplay 
of  Tumor  Vascular  Oxygenation  and  p02  in  Tumors  Using  NIRS  and  Needle  Electrode,” 
Proa  SPIE-lnt.  Soc.  Opt  Eng ,  4250,  429-436  (2001). 

Degrees  obtained  that  are  supported  by  this  award: 

Yulin  Song  was  a  Ph.D.  candidate  in  our  Biomedical  Engineering  Program  from  1995  to 
2001.  He  completed  his  Ph.D.  work  in  August  of  2001  with  a  dissertation  title  of  “Tumor 
Oxygenation  and  Physiology  in  Response  to  Therapeutic  Interventions  Investigated  by  19F  MRI 
an  NIR  Spectroscopy.”  The  current  DOD  Breast  Cancer  Grant  offered  necessary  equipment 
resources  and  intellectual  support  for  Yulin  Song  to  investigate  tumor  vascular  oxygenation 
using  Near  Infrared  spectroscopy,  while  Mr.  Song  was  also  funded  by  the  DOD  Breast  Cancer 
Research  Pre-Doctoral  Fellowship  from  1997  to  2001  (DAMD 17-97- 1-7261). 

Right  now,  Dr.  Song  is  a  post-doctoral  fellow  in  the  Department  of  Radiology  of  Stanford 
University,  California. 

Conclusions 

From  the  work  that  we  have  conducted  up-to-date,  we  can  draw  the  following  conclusions: 

1)  we  have  measured  relative  [HbOi]  changes  in  tumor  vasculature  and  tumor  tissue  p02 
under  carbogen  intervention  using  NIRS  and  needle  type  p02  electrode,  and  the  p02  data 
were  also  supported  by  the  l9F  MR  p02  mapping.  The  NIRS  data  showed  significant 
changes  in  vascular  oxygenation  accompanying  respiratory  interventions,  and  changes  in 
tumor  vascular  oxygenation  preceded  tumor  tissue  pC>2.  The  correlation  between 
normalized  tumor  A[Hb02]  and  tumor  tissue  p02  shows  that  there  is  a  threshold  of 
normalized  A[Hb02]  for  tumor  cells  to  be  efficiently  oxygenated.  This  indicates  the 
difficulty  of  oxygen  delivery  into  the  hypoxic  regions  of  tumors. 
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2)  A  combined  algorithm  can  be  developed  for  quantifying  absolute  Hb,  HbC>2,  and  Hbtotai., 
and  this  algorithm  is  particularly  made  for  our  single  source-detector  separation,  two 
wavelength,  NIR  system.  With  this  algorithm,  we  can  obtain  relatively  accurate  results  in 
[Hbtotai]  with  possible  scattering  variations. 

3)  Also,  we  learned  that  variations  in  extinction  coefficients  can  cause  significant  errors. 

References 

1 .  Mark  Cope,  “The  Application  of  Near  Infrared  Spectroscopy  to  Non  Invasive  Monitoring 
of  Cerebral  Oxygenation  in  the  Newborn  Infant”,  Doctoral  Thesis  of  Dept,  of  Medical 
Physics  and  Bioengineering,  University  College  London,  England,  (1991). 

2.  W.  G.  Zijlstra,  A.  Buursma,  H.  E.  Falke  and  J.  F.  Catsburg,  “Spectrophotometry  of 
Hemoglobin:  Absortion  Spectra  of  Rat  Oxyhemoglobin,  Deoxyhemoglobin, 
Carboxyhemoglobin,  and  Methemoglobin”  Comp.  Biochem.  Physiol.  107B,  pp.  161-166, 
(1994). 

3.  Matthias  Kohl,  Russell  Watson,  Mark  Cope.  “Optical  properties  of  highly  scattering 
media  determined  from  changes  in  attenuation,  phase  and  modulation  depth”  Optical 
Tomography  and  Spectroscopy  of  Tissue:  Theory,  Instrumentation,  Model,  and  Human 
Studies  II,  SPIE  Vol.  2979,  365-374  (1997). 

4.  Matthias  Kohl,  Russell  Watson,  Gabriel  Chow,  Idris  Roberts,  David  Delpy,  Mark  Cope. 
“Monitoring  of  cerebral  hemodynamics  during  open-heart  surgery  in  children  using  near 
infrared  intensity  modulated  spectroscopy”  Optical  Tomography  and  Spectroscopy  of 
Tissue:  Theory,  Instrumentation,  Model,  and  Human  Studies  II,  SPIE  Vol.  2979,  408-416 
(1997). 

5.  Brian  W.  Pogue,  “Error  assessment  of  a  wavelength  tunable  frequency  domain  system  for 
noninvasive  tissue  spectroscopy”,  Journal  of  Biomedical  Optics  (2000). 


List  of  Appendixes 

a)  Manuscript  submitted  to  J.  of  Biomedical  Optics. 

b)  Table  of  Contents  of  Dr.  Song’s  Dissertation. 

Among  9  chapters  in  his  dissertation,  Chapters  2,  4,  5,  6,  8  described  the  NIR  techniques 
and  methods  for  breast  tumor  measurements,  showed  results  and  drew  conclusions. 


13 


1  Mark  Cope,  “The  Application  of  Near  Infrared  Spectroscopy  to  Non  Invasive  Monitoring  of 
Cerebral  Oxygenation  in  the  Newborn  Infant”,  Doctoral  Thesis  of  Dept,  of  Medical  Physics  and 
Bioengineering,  University  College  London,  England,  (1991). 

2  W.  G.  Zijlstra,  A.  Buursma,  H.  E.  Falke  and  J.  F.  Catsburg,  “Spectrophotometry  of 

Hemoglobin:  Absortion  Spectra  of  Rat  Oxyhemoglobin,  Deoxyhemoglobin, 

Carboxyhemoglobin,  and  Methemoglobin”  Comp.  Biochem.  Physiol.  107B,  pp.  161-166,  (1994) 

3  Matthias  Kohl,  Russell  Watson,  Mark  Cope.  “Optical  properties  of  highly  scattering  media 
determined  from  changes  in  attenuation,  phase  and  modulation  depth”  Optical  Tomography  and 
Spectroscopy  of  Tissue:  Theory,  Instrumentation,  Model,  and  Human  Studies  II,  SPIE  Vol.  2979, 
365-374(1997). 

4  Matthias  Kohl,  Russell  Watson,  Gabriel  Chow,  Idris  Roberts,  David  Delpy,  Mark  Cope. 
“Monitoring  of  cerebral  hemodynamics  during  open-heart  surgery  in  children  using  near  infrared 
intensity  modulated  spectroscopy”  Optical  Tomography  and  Spectroscopy  of  Tissue:  Theory, 
Instrumentation,  Model,  and  Human  Studies  II,  SPIE  Vol.  2979,  408-416  (1997). 

5  Brian  W.  Pogue  “Error  assessment  of  a  wavelength  tunable  frequency  domain  system  for 
noninvasive  tissue  spectroscopy”, 


14 


Interplay  of  Tumor  Vascular  Oxygenation  and  Tumor  p02  Observed  Using  NIRS, 
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ABSTRACT 

Breathing  carbogen  (95%  02,  5%  C02)  to  enhance  tumor  oxygenation  may  improve  tumor  response  to 
therapy.  However,  few  studies  have  been  reported  on  the  relationship  between  tumor  oxygen  tension  (p02)  and 
blood  oxygenation  in  tumor  vasculature.  This  study  investigates  the  correlation  of  tumor  blood  oxygenation  and 
tumor  p02  with  respect  to  carbogen  inhalation  and  compares  such  correlation  with  that  in  blood.  We  used  near- 
infrared  spectroscopy  (NIRS)  to  measure  changes  of  oxygenated  hemoglobin  concentration  (A[Hb02])  and  used 
an  oxygen  needle  electrode  and  19F  MRI  for  p02  measurements  in  tumors.  The  measurements  were  taken  from 
Dunning  prostate  R3327  tumors  implanted  in  rats,  while  the  rats  were  anesthetized  breathing  air  or  carbogen.  In 
the  tissue-phantom  study,  1%  intralipid  solution  and  rabbit  blood  were  used.  NIRS  results  showed  significant 
changes  in  tumor  vascular  oxygenation  in  response  to  carbogen  inhalation.  Furthermore,  the  results  show  an 
apparent  threshold  for  tumor  p02  response.  The  slope  of  p02  versus  normalized  A[Hb02]  below  the  threshold 
obtained  from  the  tissue  phantom  was  ten  times  larger  than  the  slope  from  the  tumor,  indicating  that  tumor  cells 
are  much  harder  to  oxygenate.  This  study  demonstrates  that  the  NIRS  technology  can  provide  an  efficient,  real¬ 
time,  non-invasive  approach  to  monitoring  tumor  physiology  and  is  complementary  to  other  techniques. 

Key  Words:  Frequency-Domain  Spectroscopy,  NIR  Spectroscopy,  19F  MRI,  Tumor  Vascular  Oxygenation,  p02 
electrode,  Oxygen,  Oximetry 
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1.  Introduction 


It  has  long  been  known  that  hypoxic  tumor  cells  are  more  resistant  to  radiation  therapy  than  well-oxygenated 
tumor  cells.1  Breathing  elevated  oxygen  (100%)  or  carbogen  (95%  02,  5%  C02)  has  been  used  during  therapy  to 
improve  oxygenation.2,3  However,  little  is  known  about  oxygen  transfer  from  the  tumor  vasculature  to  tumor 
tissue  cells.  In  this  study,  we  investigated  the  relationship  between  vascular  oxygenation  and  tissue  oxygenation 
during  respiratory  challenge  in  rat  prostate  tumors  using  Near-Infrared  Spectroscopy  (NIRS),  p02  needle 
electrodes,  and  l9F  MRI  for  simultaneous  or  sequential  measurements. 

While  various  methods  are  available  for  measuring  tissue  oxygen  tension4,  many  of  them  are  highly  invasive 
precluding  the  detection  of  dynamic  changes.  Oxygen  microelectrodes5,  fiber  optic  sensors6,  and  electron  spin 
resonance7  do  facilitate  measurement  of  dynamic  changes.  MRI  has  the  further  advantage  of  providing  dynamic 
maps  of  p02,  which  can  reveal  tumor  heterogeneity.8  While  NIRS  does  not  quantify  p02,  it  can  indicate  dynamic 
changes  in  vascular  oxygenation  and  has  the  advantage  of  being  entirely  non-invasive,  providing  real-time 
measurements,  and  being  cost  effective  and  portable. 

Oxygenated  and  deoxygenated  hemoglobin  molecules  are  major  chromophores  in  tissue  in  the  near  infrared 
region  (700-900  nm),  and  they  exhibit  distinct  absorption  characteristics.  In  principle,  by  measuring  light 
absorption  and  scattering  in  tissue,  the  concentrations  of  oxygenated  hemoglobin,  [Hb02],  deoxygenated 
hemoglobin,  [Hb],  and  oxygen  saturation  of  hemoglobin,  S02,  can  be  determined  based  on  diffusion  theory. 
However,  the  theory  works  well  only  for  large  and  homogeneous  media.9, 10  Accurate  quantification  of  tumor 
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oxygenation  in  our  approach  is  currently  limited  to  relative  changes  in  [Hb02]  and  [Hb]totai  due  to  considerable 


heterogeneity  and  finite  size  of  tumors. 

The  goal  of  this  study  was  to  investigate  the  correlation  of  tumor  blood  oxygenation  and  tumor  p02  in 
response  to  carbogen  intervention  and  to  compare  such  correlation  with  that  from  regular  tissues  using  phantoms. 

2.  Materials  and  Methods 

2.1  Tumor  Model  and  p02  needle  electrode  measurements 

Dunning  prostate  rat  tumors  (three  R3327-HI  and  two  R3327-AT1)11  were  implanted  in  pedicles  on  the 
foreback  of  adult  male  Copenhagen  rats,  as  described  in  detail  previously. 12  Once  the  tumors  reached 
approximately  2  cm  in  diameter,  the  rats  were  anesthetized  with  0.2  ml  ketamine  hydrochloride  (100  mg/mL; 
Aveco,  Fort  Dodge,  1A)  and  maintained  under  general  gaseous  anesthesia  with  isoflurane  in  air  (1.3%  isoflurane 
at  1  dmVmin  air)  through  a  mask  placed  over  the  mouth  and  nose.  Tumors  were  shaved  to  improve  optical 
contact  for  transmitting  light.  Body  temperature  was  maintained  by  a  warm  water  blanket  and  was  monitored  by 
a  rectally  inserted  thermal  probe  connected  to  a  digital  thermometer  (Digi-Sense,  model  91100-50,  Cole-Parmer 
Instrument  Company,  Vernon  Hills,  IL).  A  pulse  oximeter  (model  8600,  Nonin,  Inc.,  Plymouth,  MN)  was  placed 
on  the  hind  foot  to  monitor  arterial  oxygenation  (Sa02). 

Figure  1  shows  the  schematic  setup  for  animal  experiments  using  NIRS  and  a  p02  needle  electrode.  A 
needle  type  oxygen  electrode  was  placed  in  the  tumor,  and  the  reference  electrode  was  placed  rectally.  The 
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electrodes  were  connected  to  a  picoammeter  (Chemical  Microsensor,  Diamond  Electro-Tech  Inc.,  Ann  Arbor,  MI) 


and  polarized  at  -  0.75  V.  Linear  two-point  calibrations  were  performed  with  air  (21%  02)  and  pure  nitrogen  (0% 
02)  saturated  saline  buffer  solutions  before  the  electrode  was  inserted  into  the  tumor  and  we  estimated  an 
instrumental  precision  of  1  mmHg.  P02  experimental  data  were  manually  recorded,  while  the  NIRS  data  were 
acquired  automatically.  Measurements  of  p02  and  NIRS  were  initiated  while  rats  breathed  air  for  ~10  minutes  to 
demonstrate  a  stable  baseline.  The  inhaled  gas  was  then  switched  to  carbogen  for  15  minutes  and  switched  back 
to  air. 

2.2  NIR  Spectroscopy 

NIR  spectroscopy  can  be  used  to  measure  hemoglobin  concentrations  and  oxygen  saturation  because  light 
absorptions  of  Hb02  and  Hb  are  different  at  the  wavelengths  selected  (758  nm  and  785  nm).  In  common  with  our 
previous  work13,  we  assumed  that  Hb02  and  Hb  are  the  only  significant  absorbing  materials  in  tumors  within  the 
NIR  range.  The  absorption  coefficients  comprise  the  extinction  coefficients  for  deoxyhemoglobin  (Snb)  and 
oxyhemoglobin  (sHb02)  multiplied  by  their  respective  concentrations  (Equations  1  and  2). 

Pa758  =  sHb758[Hb]  +  £Hb02  758[Hb02],  (1) 

Ea785  =  eHb785[Hb]  +  £Hb02  785[Hb02]  .  (2) 

An  NIR  homodyne  frequency-domain  photon  migration  system  allows  one  to  obtain  absorption  and  scattering 
properties  of  the  measured  sample14.  Figure  1  indicates  the  working  principle  of  the  NIR  system,  and  detailed 
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information  on  the  instrument  is  given  elsewhere13.  We  have  not  yet  completed  a  suitable  algorithm  to  compute 
pa  and  |is’  of  rat  tumors  due  to  their  finite  size  and  high  heterogeneity.  The  data  presented  here  were  analyzed 
using  only  amplitude  values  to  find  changes  in  [Hb02]  and  [Hb]totaii>  based  on  modified  Beer-Lambert’s  law 
(Equation  3).  By  manipulating  Equations  1-3,  changes  of  [Hb02]  and  [Hb],ota]  due  to  respiratory  intervention 
were  calculated  from  the  transmitted  amplitude  of  the  light  through  the  tumor  (Equations  4  and  5). 

Apa  =  Pal  -  Pali  =  LN(AB/AT)  /  L,  (3) 

A[Hb02]  =  [-5.73*LN(Ab/At)758  +  9.67*LN(AB/AT)785]  /  L,  (4) 

758  785 

A [H b>] totai  =  A[Flb]  +  A[Hb02]=[-l  ,32*LN(Ab/At)  +  6.23*LN(AB/AT)  ]  /  L,  (5) 

where  AB  =  baseline  amplitude;  AT  =  transient  amplitude;  L  =  optical  pathlength  between  source/detector;  the 
constants  given  in  the  equations  were  computed  with  the  extinction  coefficients  for  oxy-  and  deoxy-hemoglobin  at 
the  two  wavelengths  used.15  In  principle,  the  units  of  A[Hb02]  and  A[Hb]totai  are  mM,  and  L  is  equal  to  the 
source— detector  separation  multiplied  by  a  Differential  Pathlength  Factor  (DPF)16.  The  DPF  has  been  studied 
intensively  for  musclesl7and  brains18  with  approximate  values  of  4-6  and  5-6,  respectively.  Little  is  known  about 
the  DPF  for  tumors,  though  a  DPF  value  of  2.5  has  been  used  by  others19.  Since  our  focus  is  on  dynamic  changes 
and  relative  values  of  tumor  [Hb02]  with  respect  to  carbogen  intervention,  we  have  taken  the  approach  of 
including  the  DPF  in  the  unit,  viz.  modifying  eq.  (4)  as  follows: 

758  785 

A[Hb02]  =  [-5.73*LN(AB/AT)  +  9.67*LN(AB/AT)  ]/  d,  (6) 

where  d  is  the  direct  source-detector  separation  in  cm,  and  the  unit  of  A[Hb02]  in  eq.  (6)  is  mM/DPF.  We  further 
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will  show  that  our  conclusion  on  the  relationship  between  tumor  A[Hb02]  and  p02  will  not  be  affected  by  the  DPF 


after  normalizing  the  experimental  A[Hb02]  data  to  its  maximal  value. 

2.3  MRI  Instrumentation  and  Procedure 

To  support  the  findings  obtained  from  the  p02  electrode  measurements  and  NIRS,  we  conducted  MRI 
experiments  using  an  Omega  CSI  4.7  T  40  cm  system  with  actively  shielded  gradients.  A  homebuilt  tunable 
'H/I9F  single  turn  solenoid  coil  was  placed  around  the  tumor.  45  pL  of  hexafluorobenzene  (F1FB;  Lancaster, 
Gainesville,  FL)  was  administered  directly  into  the  tumor  using  a  Hamilton  syringe  (Reno,  NV)  with  a  custom- 
made  fine  sharp  (32  gauge)  needle  with  deliberate  dispersion  along  several  tracks  to  interrogate  both  central  and 
peripheral  tumor  regions  as  described  in  detail  previous6.  HFB  is  ideal  for  imaging  p02,  because  it  has  a  single 
resonance,  and  its  relaxation  rate  varies  linearly  with  oxygen  concentration.  'H  images  were  acquired  for 
anatomical  reference  using  a  traditional  3D  spin-echo  pulse  sequence.  Conventional  l9F  MR  images  were  then 
taken  to  show  the  3D  distribution  of  the  HFB  in  the  tumor.  19F  MR  images  were  directly  overlaid  over  'H  images 
to  show  the  position  of  the  HFB  in  that  slice. 

Tumor  oxygenation  was  assessed  using  FREDOM  (Fluorocarbon  Relaxometry  using  Echo  planar  imaging 
for  Dynamic  Oxygen  Mapping)  based  on  l9F  pulse  burst  saturation  recovery  (PBSR)  echo  planar  imaging  (EPI)  of 
HFB.  The  PBSR  preparation  pulse  sequence  consists  of  a  series  of  20  non-spatially  selective  saturating  90° 
pulses  with  20  ms  spacing  to  saturate  the  19F  nuclei.  Following  a  variable  delay  time  r,  a  single  spin-echo  EPI 
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sequence  with  blipped  phase  encoding  was  applied.20  Fourteen  32  x  32  PBSR-EPI  images,  with  r  ranging  from 

200  ms  to  90  sec  and  a  field  of  view  (FOV)  of  40  x  40  mm,  were  acquired  in  eight  minutes  using  the  ARDVARC 

(Alternated  Relaxation  Delays  with  Variable  Acquisitions  to  Reduce  Clearance  effects)  acquisition  protocol21.  An 

/?1(=1/T1)  map  was  obtained  by  fitting  signal  intensity  of  each  voxel  of  the  fourteen  images  to  a  three  parameter 

relaxation  model  by  Levenberg-Marquardt  least  squares  algorithm  (equation  6): 

yn  O',  j)  =  A(i,  j)-[\-(\  +  W)-  Qxp(-R\(i,  j)  •  Tn )]  (6) 

(*  =  1,2, . ,14) 

(i,j  =  1,2, . ,32) 

where  yn(i,j)  is  the  measured  signal  intensity  corresponding  to  delay  time  zn  (the  «th  images)  for  voxel  (/,/), 
A(i,j)  is  the  fully  relaxed  signal  intensity  amplitude  of  voxel  W  is  a  dimensionless  scaling  factor  allowing 
for  imperfect  signal  conversion,  and  R\(i,j)  is  the  relaxation  rate  of  voxel  (/.  /)  in  unit  of  sec'1 .  A,  W  and  R 1  are 
the  three  fit  parameters  for  each  of  the  32  x  32  voxels.  Finally,  the  p02  maps  were  generated  by  applying  the 
calibration  curve:  p02(mmHg)  =  [^(s-1)  -  0.0835]/0.001876  at  37  °C  to  the  7?1  maps.21 

2.4  Tissue  Phantom  Solution  Model 

The  tissue-simulating  phantom  comprised  a  solution  of  1%  Intralipid  in  saline  buffer.  The  oxygen  needle 
electrode,  a  pH  electrode  and  a  thermocouple  probe  (model  2001,  Sentron,  Inc.,  Gig  Harbor,  WA)  were  placed  in 
the  solution,  and  a  gas  tube  was  placed  opposite  the  NIRS  probes  to  minimize  any  liquid  movement  effects 
(Figure  2).  Source  and  detector  probes  for  the  NIRS  were  placed  in  reflection  geometry.  The  solution  was  stirred 
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at  ~  37  °C  to  maintain  homogeneity.  Fresh  whole  rabbit  blood  (2  ml)  was  added  to  the  200  mL  solution  before 


baseline  measurement.  Nitrogen  gas  and  air  were  used  to  deoxygenate  and  oxygenate  the  solution,  respectively. 

3.  Results 

3.1  Tumor  Study  Results 

Figure  3a  shows  the  temporal  profiles  of  A[Hb02]  and  p02  in  a  small  Dunning  prostate  R3327-H1  tumor  (1 .5 
cm2)  with  respiratory  challenge  measured  simultaneously  with  N1RS  and  the  p02  needle  electrode.  After  a  switch 
from  air  to  carbogen,  A[Hb02]  increased  rapidly,  along  with  tumor  tissue  p02.  Figure  3b  was  obtained  from  a 
second  large  prostate  tumor  (3.1  cm2)  and  now  the  NIRS  response  was  biphasic,  while  the  electrode  showed  a 
slower  p02  response.  In  a  third  tumor  NIRS  behaved  as  before,  but  p02  did  not  change  (Fig.  3c). 

In  two  tumors  from  a  separate  subline  (Dunning  prostate  R3327-AT1),  NIRS  and  l9F  MRI  were  taken  with 
carbogen  challenge  sequentially  (Fig.  4).  NIRS  response  showed  vascular  oxygenation  changes  as  before,  and 
FREDOM  revealed  the  distinct  heterogeneity  of  the  tumor  tissue  response.  Initial  p02  was  in  the  range  of  1  to  75 
mmHg  and  carbogen  challenge  produced  values  in  the  range  6  to  360  mmHg.  Representative  voxels  are  shown. 

3.2  Tissue  Phantom  Solution  Study  Results 

Figure  5  shows  a  similar  temporal  profile  for  A[Hb02]  and  p02  measured  from  the  tissue  phantom  during  a 
cycle  of  gas  change  from  air  to  nitrogen  and  back.  The  first  five  minutes  were  measured  as  a  baseline  after 
adding  2  ml  blood.  Bubbling  nitrogen  deoxygenated  the  solution  and  caused  the  p02  values  to  fall;  A[Hb02] 
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declined  accordingly  with  a  small  time  lag.  After  the  bubbling  gas  was  switched  from  nitrogen  to  air,  both 


A[Hb02]  and  p02  started  to  increase  simultaneously,  but  the  recovery  time  of  A[Hb02]  to  the  baseline  was  faster 
than  that  of  p02.  The  small  time  lag  between  the  changes  of  A[Hb02]  and  p02  is  probably  due  to  the  allosteric 
interactions  between  hemoglobin  and  oxygen  molecules.  According  to  the  hemoglobin  oxygen-dissociation 
curve22, 2\  oxyhemoglobin  starts  to  lose  oxygen  significantly  when  p02  falls  below  70  mmHg  at  standard 
conditions  (pH=7.4  and  temperature=37  °C).  The  same  principle  can  explain  why  A[Hb02]  has  a  faster  recovery 
than  that  of  p02.  Namely,  A[Hb02]  became  saturated  as  p02  increased  from  0  to  70  mmHg,  while  the  solution 
was  being  oxygenated.  Importantly,  A[Hb]tota|  remained  unchanged,  as  expected,  during  the  cycle  of 
deoxygenation  and  oxygenation. 

3.3  Correlation  of  p02  and  A[Hb02]/A[Hb02]max 

To  compare  the  data  from  tumors  and  the  phantom  for  the  relationship  between  changes  in  tumor  vascular 
oxygenation  and  tumor  tissue  p02,  we  normalized  A[Hb02]  at  the  maximum  value  during  the  carbogen 
intervention,  i.e.,  A[Hb02]/A[Hb02]max.  One  advantage  of  using  normalized  A[Hb02]  is  to  eliminate  the 
uncertainty  of  DPF. 

Figure  6a  replots  the  data  given  in  figures  3a  and  3b,  showing  a  direct  relationship  between  the  normalized 
A[Hb02]  and  p02  in  the  tumors  measured  with  NIRS  and  the  p02  needle  electrode.  It  shows  that  p02  in  the 
tumors  varied  linearly  with  a  slope  of  4. 1  in  case  1  and  a  slope  of  3.6  in  case  2  before  A[Hb02]  reached  80  %  and 
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60  %  of  A[Hb02]max  in  cases  1  and  2,  respectively.  Thereafter,  tumor  p02  increased  sharply,  indicating  a  threshold 


for  tumor  tissue  p02  to  increase.  While  NIRS  results  tended  to  be  similar  for  several  tumors,  p02  electrode 
measurements  showed  considerable  variation  even  in  the  same  tumor  type,  suggesting  distinct  tumor 
heterogeneity.  This  was  substantiated  by  the  l9F  MR  p02  mappings:  indeed,  in  some  cases,  p02  values  did  not 
change  with  respiratory  challenge,  especially  when  baseline  p02  values  were  lower  than  10  mmHg. 

Figure  6b  shows  the  relationship  between  normalized  A[Hb02]  and  p02  measured  from  the  phantom.  The 
threshold  occurred  at  essentially  the  same  value  of  A[Hb02]/A[Hb02]max  (-80%)  as  observed  for  tumor  case  2. 
However,  the  phantom  showed  a  much  larger  slope  between  p02  and  A[Hb02]/A[Hb02]max,  i.e.,  39.5,  prior  to  the 
threshold. 

Figure  6c  summarizes  the  data  obtained  from  the  tumors  and  phantom  and  plots  them  with  reversed  axes.  In 
the  phantom,  p02  had  increased  from  0  mmHg  to  70  mmHg  by  the  stage  that  A[Hb02]  reached  to  90%  of 
A[Hb02]max.  This  behavior  is  reminiscent  of  the  hemoglobin  oxygen  dissociation  curve.  By  comparison,  the 
tumor  p02  increased  only  from  14  mmHg  to  30  mmHg  in  case  1  and  from  10  mmHg  to  21  mmHg  in  case  2  when 
A[Hb02]  approached  to  90  %  of  A[Hb02]max. 

4.  Discussion  and  conclusion 

Tumor  oxygenation  involves  a  complex  interplay  of  multiple  compartments  and  parameters:  blood  flow, 
blood  volume,  blood  vessel  structure,  and  oxygen  consumption.  NIRS  provides  a  global  non-invasive  estimate  of 
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average  vascular  oxygenation  encompassing  arterial,  venous  and  capillary  compartments.  In  agreement  with  our 
previous  observations13,  the  A[Hb02]  response  is  often  biphasic,  which  we  believe  represents  rapid  elevation  of 
arterial  oxygenation,  followed  by  more  sluggish  capillary  components.  Comparison  with  electrode  measurements 
indeed  revealed  that  tumors  are  heterogeneous.  Like  an  NIRS  measurements,  p02  electrodes  provide  rapid 
assessment  of  p02  facilitating  real  time  observation  of  dynamic  changes.  In  Fig.  3a,  p02  starts  at  a  baseline  value 
~  15  mmHg  and  increases  rapidly  in  response  to  respiratory  challenge  with  carbogen.  Indeed,  the  rate  approaches 
that  of  the  vascular  compartments.  In  a  second  tumor,  where  the  interrogated  location  showed  a  slightly  lower 
p02,  the  tissue  response  was  more  sluggish  (Fig.  3b).  For  a  third  HI  tumor,  local  baseline  p02  was  found  to  be  <  5 
mmHg  and  this  did  not  changed  with  carbogen  inhalation  despite  the  response  observed  by  NIRS.  This  suggests  a 
danger  of  comparing  a  global  vascular  measurement  with  regional  tumor  p02,  since  tumors  are  known  to  be 
highly  heterogeneous. 

Indeed,  FREDOM  measurements  in  Fig.  4  revealed  the  heterogeneity  in  baseline  oxygenation  within 
individual  tumors  of  this  second  tumor  subline.  Baseline  p02  ranged  from  1  -  75  mmHg  and  response  to 
carbogen  was  variable  in  terms  of  rate  and  extent,  as  seen  for  the  HI  subline  using  electrodes.  As  with  the 
electrodes,  the  better  oxygenated  tumor  regions  showed  a  faster  and  greater  response  to  carbogen  inhalation. 

The  phantom  measurements  indicate  the  reliability  of  the  NIRS  technique.  Both  phantom  and  tumor 
measurements  indicate  a  sharp  increase  of  p02  after  the  normalized  A[Hb02]  reaches  about  80%  of  its  maximum. 
This  coincides  with  the  recognized  sigmoidal  binding  of  oxygen  to  hemoglobin.  The  difference  between  the 
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phantom  and  tumor  results  reveals  that  the  p02  response  to  [Hb02]  changes  in  tumor  is  about  10  times  smaller 
than  that  in  tissue  phantom  when  p02  values  are  relatively  low  (0  to  30  mmHg).  This  demonstrates  a  greater 
hurdle  to  delivering  oxygen  from  blood  vasculature  into  the  hypoxic  regions  of  the  tumor  than  into  regular  tissues. 
The  slope  at  low  A[Hb02]/A[Hb02]max  in  the  tumors  presumably  also  reflects  oxygen  consumption  by  the  tumor 
tissues,  as  oxygen  initially  becomes  available. 

Both  NIRS  and  electrodes  offer  essentially  real  time  measurement  of  changes  in  oxygenation,  which  can  be 
rapid  (Fig.  3).  Indeed,  the  inflow  kinetics  of  vascular  02  detected  by  NIRS  are  similar  to  those  previously 
reported  in  the  HI  tumor  line  following  a  bolus  of  the  paramagnetic  contrast  agents  Gd-DTPA24.  FREDOM  has 
lower  temporal  resolution,  but  reveals  the  tumor  heterogeneity  and  differential  response  of  regions  exhibiting 
diverse  baseline  p02.  The  results  here  correspond  closely  with  more  extensive  observation6,  8’  2I.  While 
FREDOM  currently  requires  6V2  minutes  per  p02  map,  we  have  previously  demonstrated  alternative  data 
acquisition  protocol  achieving  1  s  time  resolution  in  a  perfused  heart,  albeit  providing  less  precision  in 
measurements  and  only  a  global  determination.25 

In  conclusion,  we  have  measured  relative  [Hb02]  changes  in  tumor  vasculature  and  tumor  tissue  p02  under 
carbogen  intervention  using  NIRS  and  needle  type  p02  electrode,  and  the  p02  data  were  also  supported  by  the  l9F 
MR  p02  mapping.  The  NIRS  data  showed  significant  changes  in  vascular  oxygenation  accompanying  respiratory 
interventions,  and  changes  in  tumor  vascular  hemoglobin  oxygenation  preceded  tumor  tissue  p02.  The  correlation 
between  normalized  tumor  A[Hb02]  and  tumor  tissue  p02  shows  that  there  is  a  threshold  of  normalized  A[Hb02] 
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for  tumor  cells  to  be  efficiently  oxygenated.  This  indicates  the  difficulty  of  oxygen  delivery  into  the  hypoxic 


regions  of  tumors.  This  study  demonstrates  that  the  NIRS  technology  can  provide  an  efficient,  real-time,  non- 
invasive  approach  to  monitoring  tumor  physiology,  and  emphasizes  the  need  to  develop  an  imaging  technique  to 
reveal  spatial  heterogeneity. 

Acknowledgement 

This  work  was  supported  in  part  by  The  Whitaker  Foundation  (HL),  NIH  ROl  CA79515  (RPM),  and  the 
Department  of  Defense  Breast  Cancer  Initiative  grant  BC990287  (HL).  NMR  experiments  were  conducted  at  the 
Mary  Nell  &  and  Ralph  B.  Rogers  MR  center,  an  NIH  BRTP  Facility  no.  5-P41-RR02584.  We  are  grateful  to 
Sherwin  Yen  and  William  McSweeny  for  assistance  with  data  acquisition  and  Dr.  Peter  Peschke  for  providing  the 
original  tumor  cells. 


13 


Figure  Captions 


Figure  1.  Schematic  experimental  setup  of  one  channel,  near  infrared,  frequency  domain  IQ  instrument  for  tumor 
investigation  in  vivo.  The  5-mm  diameter  fiber  bundles  deliver  the  laser  light,  comprising  two  wavelengths  (758 
and  785  nm)  and  detect  the  laser  light  following  absorption  through  the  implanted  tumor.  The  p02  needle 
electrode  measures  tumor  tissue  p02. 

Figure  2,  Experimental  setup  for  phantom  study  using  1%  Intralipid  in  saline  buffer.  NIRS  probes  were  placed  in 
reflectance  mode,  while  the  gas  bubbler  was  placed  opposite  to  minimize  liquid  movement  effects.  After  adding 
2  ml  of  rabbit  blood  to  a  200  ml  solution,  nitrogen  gas  and  air  were  introduced  to  deoxygenate  and  oxygenate  the 
solution,  respectively. 

Figure  3.  Simultaneous  dynamic  changes  of  A[Hb02]  and  p02  in  R3327-HI  rat  prostate  tumors  using  NIRS  and 
p02  needle  electrode;  (a)  A  smaller  tumor  (1.5  cm3)  showed  a  rapid  p02  response  (Case  1),  whereas  (b)  a  bigger 
tumor  (3.1  cm3)  showed  a  slower  p02  response  (Case  2).  (c)  In  a  third  tumor  where  regional  baseline  p02  was  <  5 
mmHg,  there  was  no  p02  response  (Case  3).  The  unit  of  A[Hb02]  is  mM/DPF,  where  DPF  is  equal  to  the  optical 
path  length  divided  by  the  source-detector  separation.  Dotted  vertical  line  marks  the  time  when  the  gas  was 
changed.  The  oxygen  needle  electrode  was  placed  in  tumors  where  p02  had  a  baseline  value  of  about  10  mmHg. 


14 


Figure  4.  Dynamic  changes  of  A[Hb02]  and  p02  in  R3327-AT1  rat  prostate  tumors  measured  using  NIRS  and  l9F 
MR  p02  mapping.  The  solid  curve  represents  A[Hb02]  and  the  solid  line  with  solid  circles  represents  mean  p02  ± 
SE  (Standard  Error)  of  21  (Fig.  4a)  and  45  (Fig.  4b)  voxels  of  each  tumor.  Dashed  lines  are  4  representative 
voxels  for  each  case.  After  a  gas  switches  from  air  to  carbogen,  both  the  mean  p02  of  each  tumor  increased. 
Individual  voxels  showed  totally  different  response,  indicating  heterogeneity.  The  tumor  sizes  were  3.2  cm3  and 
2.7  cm3  for  (a)  and  (b),  respectively. 

Figure  5.  Simultaneous  dynamic  changes  of  A[Hb02],  A[Hb]t0,ai  and  p02  in  the  solution  phantom  measured  using 
NIRS  and  p02  needle  electrode.  The  dark  solid  curve  is  for  A[Hb02],  the  lighter  solid  line  is  for  A[Hb]totai,  and  the 
solid  circles  show  p02  values  in  the  phantom  solution.  After  ~3  minutes  of  baseline,  the  bubbling  gas  was 
changed  from  air  to  nitrogen  to  deoxygenate  the  solution  and  then  switched  back  to  air  to  reoxygenate  the  solution. 
The  unit  of  A[Hb02]  is  mM/DPF. 

Figure  6.  Changes  of  tumor  tissue  p02  with  normalized  changes  of  oxygenated  hemoglobin  (a)  in  tumors 
measured  with  NIRS  and  a  p02  needle  electrode,  and  (b)  in  the  phantom  solution  using  the  same  NIRS  and  p02 
needle  electrode.  The  slope  between  p02  and  A[Hb02]/A[Hb02]max  is  3.6  when  A[Hb02]/A[Hb02]max  is  less  than 
0.8  in  Case  2  tumor,  while  the  slope  obtained  from  the  phantom  solution  is  39.5,  which  is  ten  times  larger  than 
that  in  Case  2.  (c)  Summarizing  of  relationship  between  p02  and  normalized  A[Hb02]  obtained  from  both  tumors 
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and  the  tissue  phantom.  Dotted  horizontal  line  shows  90%  of  normalized  A[HbC>2]. 
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